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c/3 ■ Abstract 

O 

In this review paper, we present and critically re-examine the formal ex- 
pressions for the electrophoretic effect and the ionic field appearing in the un- 
published Yale University PhD dissertation of W. S. Wilson which form the 
basis of the Onsager- Wilson theory of the Wien effect in the binary strong elec- 
trolyte solutions. It is pointed out that some of the integrals that make up the 
flow velocity formula obtained in the thesis and he evaluated at the position 
of the center ion in the ionic atmosphere (i.e., the coordinate origin) diverge. 
f***3 \ Therefore they cannot be evaluated by means of contour integrals in the man- 

C*~) ■ ner performed in his thesis for the reason pointed out in the text of this paper. 

**) | In this paper, the results for the integrals in question are presented, which are 

l/~) ■ alternatively and exactly evaluated. The details will be described in the follow- 

up paper presented elsewhere together with the improved formula for the Wien 
effect on conductivity. 



I. Introduction 

Non-Ohmic phenomena in charged media, such as plasmasP semiconduc- 
tors,^ ionic solutions, 3 subjected to high external electromagnetic fields are much 
observed in recent years. They not only have numerous practical applications 
in science and engineering, but also pose much challenge theoretically. Espe- 
cially, in small systems^ confined in space, which consequently are subject to 
a very large field gradient, nonlinear phenomena are usually unavoidable and 
their theoretical treatment has been of considerable interest as a large body of 
literature on the subject matter attests to the interest, relevance, and activity in 
the fields. The usual linear theory would not suffice for an adequate theoretical 
treatment, and one would require some sort of nonlinear theory of irreversible 
processes, whether a kinetic theory (statistical mechanics approach) or macro- 
scopic irreversible thermodynamics approach is taken. The present author has 
developed a nonlinear theory of irreversible processes both in the kinetic theory 
(molecular theory) approach^ or the phenomenological irreversible thermody- 
namics approach, 7 when the medium (liquids or gases) is electrically neutral. 
Naturally, extension of the theories mentioned of neutral media has been one of 
the aims of recent research of the present author. 

In electro-physical chemistry, non-Ohmic behavior of conductivity was ob- 
served since 1920s. It is known as the Wien effectPin electrolyte solutions. The 
theories of electrical conductivity in electrolyte solutions have been developed 
by Lars Onsager and his collaborator d 9 * 10 * 11 ! 12 ! over a period stretching over 
twenty some years beginning from the late 1920s. Onsager's theories elucidate 
the important effects underlying the conduction phenomena in charged media 
which we may make use of for many conduction phenomena observed in the 
fields and matters mentioned of current interest. Despite their significance and 
implications there has been very few follow-up studies of his theories by other 



authors in which the theories could have been critically reviewed and analyzed 
for further development at the basic level. Although, superficially, they seem 
dated from the current standpoint it appears there still are many features and 
techniques in them from which we could much benefit for the current new things 
we daily encounter in the fields mentioned earlier. On the strength of this line 
of reasoning and motivated by the possibility of utility of his theories for new 
related problems, the present author has been studying his theories as closely as 
possible at the basic level. We would like to report that even an improvement of 
one of the results on the Wien effect in binary electrolytes can be achieved that 
would affect the experimental interpretations in the past of the conductivity in 
electrolyte solutions carried out on the basis of the Wien effect on electrolytes. 
In the path breaking paper of ionic conductivity in strong electrolyte solu- 
tions subject to an external electric field, Lars Onsagep^l in 1927 derived a linear 
conductance formula for strong electrolyte solutions. To derive the formula he 
took into consideration the electrophoretic effect and the ionic held relaxation 
effect, which were estimated by means of physical arguments — on the basis of 
the Stokes law and a relaxation time argument. The formula derived estab- 
lished the limiting laws of ionic conduction, which ever since has appeared in 
textbooks of physical chemistry and, particularly, physical chemistry of elec- 
trolyte solutions. This early linear theory was followed up by W. S. Wilson^ in 
his 1936 Yale PhD thesis under Onsager's supervision in which he obtained the 
formula for ionic conductance in binary strong electrolytes in high external elec- 
tric fields — the Wien effect, which was experimentally discovered by M. WierP 
in 1920s. His thesis has not been published in a journal, but a summary and the 
tables for the functions representing the electrophoresis effect and the ionic field 
relaxation time effect have appeared in the well-known monumental reference 
book on physical chemistry of electrolyte solutions by Harned and OwenP^I The 



reason for not publishing Wilson's results, according to Harned and Owen, was 
the partial Fourier transform method used in Wilson's thesis should be replaced 
by a full Fourier transform method to be developed later. However, it took two 
decades for such a method to appear in 1957. In the intervening period, the 
aforementioned tables prepared by Wilson have been used in connection with 
the Wien effect in the literature. 

When compared with experiment, the Wien effect formula obtained, how- 
ever, does not fully account for the observed conductance data on strong elec- 
trolytes. Therefore, the theory should be improved, but an improvement cannot 
be achieved unless we know the theory in its intimate detail. Because of the 
aforementioned deficiency, the ionic conductance formula for the Wien effect has 
been seminal and played an important role in the development of the idea of 
the ionic association theory put forward by A. Patterson and his collaborators^ 
who have made use of Bjerrum's ionic association theory 16; which, coupled with 
the weak electrolyte theory of OnsagerpS has enabled them to account for the 
observed conductance dataP^l Thus the idea of degree of ionic association took 
a place in physical chemistry lore. 

Having formulated a nonlinear transport theory — generalized hydrodynamic 
theory — of non-electrolytes (neutral liquids and gases) by means of kinetic the- 
ory^ and its phenomenological version^ the present author in the recent years 
has been investigating application of the aforementioned generalized hydrody- 
namic theory to electrolytic conduction under high external electric fields often 
encountered in various small systems^ of recent interests in various fields. The 
natural course to take therefore has been to learn the Onsager line of theo- 
ries before deciding to apply them to phenomena in mind. Unsatisfied by the 
brief summary of the theory of Wien effect given by Wilson in the monograph of 
Harned and Owen, 14 the present author decided to have a close look at Wilson's 



thesis. 

In his thesis, by solving the governing equations for the distribution functions 
and Poisson equations which were earlier established by Onsager and Fuossp^H 
Wilson obtains the solutions for the distribution functions and ionic potentials. 
Therewith he then obtains a formal expression for the ionic field. He also solves 
the Navier-Stokes equation for velocity by a (partial) Fourier transform method 
and obtains a formal expression for the velocity field. These formal expressions 
(solutions) are made up of complicated integrals [in fact, cos(ax) and sin(ax) 
transforms, where x is the axial coordinate and a the transform variable] in- 
volving zeroth-order Bessel functions of second kind Ko(Xp) with non-simple 
argument where p is the radial distance from the central ion in the cylindri- 
cal coordinate system employed and A is an irrational function of a, the field 
strength, and the Debye length characterizing the ion atmosphere of the elec- 
trolyte solution of interest. Since the physical quantities of interest are those 
of the central ion located at the origin of the coordinates (the center of ionic 
atmosphere), he elects to take x = and p = within the integrals before 
evaluating them. Then he resorts to contour integration methods and manages 
to evaluate the integrals in closed analytic forms. The results give rise to the 
well-known formula for the Wien effect, which was later found to give qualita- 
tively correct results compared with experiments, but not quantitatively in the 
high field regime. 

The aforementioned step of taking special values of the variables, x = 
and p = 0, is found to give rise to a difficulty because a couple of the resulting 
integrals preclude the contour integral method that involves a contour at infinite 
radius in the complex plane. In fact, when closely examined, they diverge as 
the argument (a) tends to infinity on the positive real axis. Therefore, it is not 
permissible to take p = 0, which allows to take the Bessel function Ko(Xp) at p = 



before evaluating the integral. This means, mathematically, that the integrals 
are not uniformly convergent with respect to p because of the characteristics of 
Ko(\p) at p — in the integrals. In the follow-up 19 of this tutorial review article, 
the integrals appearing in the formal expressions for the velocity field and ionic 
field in Wilson's thesis are evaluated alternatively and exactly. The integrals 
can be reduced to definite integrals of another zeroth-order Bessel function of 
second kind, Iq (Xp), along the imaginary axis, which may be evaluated term 
by term in terms of quadratures of algebraic functions. They are rather easily 
amenable to numerical integration. The main purpose of this tutorial article is 
to introduce to the reader the solution procedures leading to Wilson's results 
from the governing equations for the distribution functions and the Navier- 
Stokes equation with the local force field constructed from the solutions of the 
governing equations. They are not simple at all, yet not readily available in the 
main stream scientific literature easily accessible at present. (Wilson's thesis 
is only available on interlibrary loan from Yale University.) Besides, Wilson's 
dissertation contains a number of important typographical errors often requiring 
checking the whole solution processes to correct them. In view of the importance 
and significance of the Wien effect in nonlinear physico-chemical processes in 
electrolytes it is hoped to provide the reader with the details of the Onsager- 
Wilson theory of electrolytes. This author believes that their theory should 
provide significant insights and methodologies for the modern investigations in 
electrical conductivities and transport properties of electrically charged fluids 
and condensed matter physics. 

This article is organized as follows. In Sec. II, in view of the fact that 
Wilson's thesis was not published and its details are not available to the general 
audience in the mainstream scientific literature it seems appropriate to provide 
the essential materials beyond what is presented in the monograph of Harned 



and Owen. As a matter of fact, it contains valuable lessons and insights for 
us to theoretically and mathematically treat physical phenomena in electrolytes 
and plasmas subjected to external electric fields. Thus, for notational purposes 
we will first present the governing equations and the formal solutions of the 
governing equations. Then the formal solution is obtained for the Navier-Stokes 
equation by using the local force field formula supplied by the aforementioned 
solutions for the governing equations. The velocity field and ionic field obtained 
are formal, being given in terms of one-dimensional Fourier transforms of zeroth- 
ordcr Bcssel functions mentioned earlier. We then will discuss the procedure 
used by Wilson to evaluate them, for which he uses the method of contour 
integration and residues. Instead of his contour integral methods we will use 
our own contours different from his. They are more easily comprehensible in our 
opinion. We will indicate an alternative method of evaluation of the integrals, 
but only the results are presented, the details of which will be dealt with in the 
regular article by the present author elsewhere. This method gives rise to finite 
integrals of the zeroth-order Bessel functions of second kind, Jo (z), along the 
imaginary axis in the a plane, which in fact can be expressed as quadratures 
of algebraic functions — polynomials. In Sec. IV, Wilson's results obtained for 
the velocity field and the ionic field are then used for the formula for ionic 
conductance of binary strong electrolytes. Sec. V is for the concluding remarks. 

II. The Governing Equations and Solutions 

To treat transport phenomena in electrolyte solutions Onsager and Fuosa^l 
proposed a Fokker-Planck type equation for ion pair correlation functions — 
probability distribution functions. The distribution functions are introduced in 
the following manner. 

Let the two volume elements dV\ and dVi be at V2 and ri, respectively, at 



time (duration) t. Let the fraction of time tj of t be the time duration in which 
particle j is found in dV\ . The probability (density) rij of finding j in dV± during 
that time duration is then, by the ergodic hypothesis p3 

^ =n 1 dV 1 , 
t 3 

whereas the probability n% of finding i in dV2 during that time duration is 

— = riidVi. 

t 

Now let tji be the time duration of tj (i.e., the fraction of tj) in which an ion i 
is found in dV2 , given an ion j in dV\ . The conditional probability riji of finding 
an i ion in dVi is then 



f- = rijidV-2. 



Similarly, in the reversed case, 



Since 



f- =n ij dV 1 . 



tji "ij 



it follows 



We may write 



njrijidVidV2 = ninijdV\dV2- (1) 



fji (ri,r 2 i) = rijUji = UiUij = fy (r 2 , r 12 ) . (2) 

Here ri and r 2 are position vectors of ion j and ion i, respectively, in a suitable 
coordinate system and 

r = r 2i = r 2 - r! = -r i2 . (3) 

The fji (ri,r 2 i) is the probability of finding ion i at r 2 in dV 2 at distance r 2 i 
from ion j at ri in dV\ , and inversely for fy (r 2 , ri 2 ) . These are pair distribution 
functions of the particle pair (i,j). 



The charge density pj at a position r from ion j is then given by 

f.. 

i i ■* 

where e, is the charge on ion i. The potential due to the j ion and its atmosphere 

is denoted tpj (ri, r 2 i), which obeys the Poisson equation 

n u( , s 4tt ^^ /jt(ri,r 2 i) , . 

V-V^-(r 1 ,r 2 i) = -— 2_s^—^ tn, (5) 

i 3 

where D is the dielectric constant. 

Let the velocity of ion i in the neighborhood of ion j be Vji — Vji (ri , r 2 i) and 
inversely v, 3 = v,y (r 2 ,ri 2 ). With the relative coordinate r introduced, the po- 
sition coordinates ri and r 2 can be suppressed and we may write Vji (ri, r 2 i) = 
v^ (r) and Vy (r 2 , r i2 ) = v^ (-r). 

In the Brownian motion modeP^ it is assumed that 

Vji (r) = v (r) + uh (Kji -k B T\n fe) , (6) 

Vij (-r) = v (r) + Uj (K t j - k B T\n f tj ) , (7) 

where v (r) is the barycentric velocity, the total force K^ is given in terms of 
the external force F, and fluctuating forces of Brownian motion 

Kji = ki + fluctuating forces, 

and w, is the mobility of ion i that produces velocity w^ki when k^ acts on 
ion i. In the relative coordinates adopted, the probability /^(ri,r 2 i) — pair 
distribution function — is assumed to obey the equation of continuity 



dfji (ri,r 2 i) 



= Vi • [fjiVji (n, r 2 i) + fijVij (r 2 , n 2 )] . (8) 



at 

With the Brownian motion model (JSJ and ([7]) a steady state Fokker-Planck type 
equation is obtained from Eq. ([5} for fjf. 

-V • [fjiVji (r) - f^v^ (-r)] = Vi ■ Wj (/«K« - *i,TVi/y) (9) 

+ V 2 • w t (fjiKji - ksTVifji) . 



Equations ([5]) and © are the basic steady-state governing equations for the 
distribution functions fji and potentials ipj in the Onsager-Fuoss theory of 
electrolytic transport processes in external electric fields. They are coupled as 
will be clear presently. Further assumptions/approximations are made to Eq. 

©• 

Since the fluctuating forces making up the total force Kjj mainly originate 
from ionic interactions within the ionic atmosphere and between ionic atmo- 
spheres of ions i and j , they must be getting contributions from ionic potentials 
ipi and ipj. Neglecting the ionic forces of O (e 2 ) or higher, the total forces in 
the Onsager-Fuoss theory K^j are looked for in the form 

Kji =ki- eiVitp'i (0) - aVi^j (ri,r 21 ) , (10) 

where the prime denotes the nonequilibrium correction to ipi beyond the Debye- 
Huckel equilibrium potential im obeying the Poisson-Boltzmann equation 

i— 1 

Here the Debye parameter (inverse Debye length) k, is defined by 



\ 



' JZ^el (12) 



k B TD . 

1=1 



Thus we now see that Eqs. ([5]) and (© are coupled. We look for the solutions 
in the forms 

^•(r 1 ,r 21 )=^°(r)+^(r 1 ,r 21 ), (13) 

riji (ri,r 2 i) = n% (r) +n' Jt (ri,r 2 i) . (14) 

On substitution of Eq. (TlU)) . (fl^j) . and (1X4^1 . the steady-state governing equations 
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for fji are now given by the equation 



V • [vji (r) - Vy-(-r) ] f^ =Uiki ■ Vfy (r) - ojjkj ■ Vfa (r) 

k B T. 



V • ujiTijni 



V ■ LOjUjUi 



eiViPj (r) + -^—^fji (r) 
ejVVi (-') + — V/,-<(r) 



(15) 



(m = 1,2,- ■• ,s) 



Note that Vj • k^ = is made use of in this equation since in the experiments 
of interest the external forces ki are maintained constant in space. For this 
equation we have also made use of the fact that 



fji - mrij = O (etej) = 0(e 2 ), 

fji ( r ) = fij (- r ) > 



(16) 
(17) 



and neglected the terms containing VV^(O) and V?/>'(0) owing to the fact that 
they are 0(e 2 ). This final form of the governing equation for fji are coupled to 
the Poisson equation 



VV,(r) = -^E 



47r -A fji (r) e H 



d ^ m 

i=\ 3 



(18) 



Wilson solved this set for binary electrolyte solutions on which the external 
electric field X acting in the direction of positive x axis: 



fci = eiX, k 2 = e 2 X. 



(19) 



Written out for a binary electrolyte — in fact, a uni-uni electrolyte, Eq. (fT5|) and 
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Poisson equations becomes a set of six coupled equations: 

n 2 e { V 2 [V>i (r) + ^ (-r)] } + 2k B Tf n (r) = 0, (20) 

-n 2 e { V 2 [V> 2 (r) + ^a (-r)] } + 2fe B T/ 22 (r) = 0, (21) 
Xe (wi + w 2 ) V a /i 2 (r) - n 2 e {wxV 2 ^ (-r) - w 2 V 2 ^i (-r)} 

+k B T (wi + w 2 ) V 2 / 12 (r) = 0, (22) 
-Xe (wa + wi) Vx/ai (r) - n 2 e {wxV 2 ^ (-r) - WaV^i (r)} 

+fc B T ( Wl + w 2 ) V 2 / 21 (r) = 0, (23) 

Afrp 

VVi(r)--[/n(r)-/ 12 (r)], (24) 

4-7TP 

VV 2 (r) = -— [/ 21 (r)-/ 22 (r)]. (25) 

Here V x = d/dx. Note that for a binary electrolyte 

n 1 ei+n 2 e 2 ^0 (26) 

by electroneutrality and n\ = n 2 . These coupled equations are solved subject 
to the boundary conditions. At this point, it will be useful to remind ourselves 
that these governing equations, when solved, would provide us with the electrical 
perturbations to the pair distribution functions and the potentials. In Onsager's 
theory there is no provision for calculation of the equilibrium part of the pair 
distribution functions in the absence of the external field. The equilibrium 
distribution functions should be calculated and supplied separately, if necessary. 
For them we may make use of the modern theory of equilibrium charged fluidsP^I 

III. Boundary Conditions 

Eqs. (|2"0"l) - (f2"5")) are subject to the following boundary conditions. 
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A. Flow 

The number of ions Fij (Q) leaving the interior of a region S is 

8F- r 

~~m = J s fji (r) {e " ■ [Vjl (r) " V?J (_r)]} dS 

= / V • fji (r) [ Vji (r) - Vij (-r)] dO = 0, (27) 

where e„ is the unit vector normal to the surface S. It is equal to zero by the 
law of mass conservation. Since volume £2 can be as small as possible, the limit 
of fl — > may be taken: 

lim f f Tl (r) e n • [v ri (r) - v« (-r)] dS = 0. (28) 

S2 ^° Js 

This means that the flow field fji (r) [vjj (r) — v^ (— r)] must be sourceless. 

That is, /jj (r) must satisfy this condition. According to Eq. ([15)) this flow 

condition implies 

fji (r) [Vji (r) - v^ (-r)] = w t [kj^ (r)] - Uj [k 3 f l3 (-r)] 

- [wieiriiiijVipj (r) + UjejmnjVipi (— r)] 
-fcaTCwi + w^V/iiCr). (29) 

This condition will be made use of later. 

B. Ionic Fields 

The space charge within £7 is 

/ dflpj (r) = — — / dftV • V^j (r) (30) 

Jo 47r ,/ n 

As fi —> 0, the volume of the ion j at the origin, the left-hand side of this 
equation is equal to ej if the charge is located at the origin of the coordinates, 
or equal to zero if the charge is not in Q. Thus we may write it as 

Pi (r) = e^(r). (31) 

13 



Therefore 



or 



where 



4-7TP ' 

V-V^(r) = — ff-5(r) (32) 



a /dS.V^(r) = -^ (33) 



(5=1 if ej is at r = 
= otherwise. 

Since near r = 0, the screening effect is equal to zero and hence the potential is 
equal to 



and in the large distances 



^ = Wr (34) 



lim V, = 0. (35) 

r— foo 



Hence the boundary conditions on ipj are deduced to be 

ipj - ^- < oo as |r| -> 0, (36) 

V>j (oo) = 0. (37) 

If the charge is of finite size Gj , then this should be modified to 

g - 
ipj - — p < oo as |r| > <jj, (38) 

ipj (oo) = 0. (39) 

The Onsager- Wilson theory does not take the case of finite ion sizes into con- 
sideration as should have been. 
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C. Symmetry Conditions for Potentials and Distribution 
Functions 

1. Potentials 

Since the distribution functions are even 

/« (r) = f u (-r) , f n (r) = f n (-r) , (40) 

fii M = fij (-r) , (41) 

it is possible to deduce the symmetry conditions for the potentials. From the 
Poisson equations we find 



4?re, 

Dn 

Aire 



V 2 ^- (r) = —=B [fa (r) -/*(!•)], (42) 



V 2 ^ (~r) = -^ [/« (- r ) - & (- r )] • (43) 



Adding the two equations yields 

,2r. - ■ - - 47re 



V 2 [^ (r) + V* (-r)] = - — fe/tf (r) + */« (r)] • (44) 



From this follows 



Aire 



V 2 [^ (-r) + ^ (r)] = ~ [Zjfjj (r) + zj u (r)] . (45) 

Here e, = ezi and ej = ezj. Subtracting it from the previous equation, we 
obtain 

V 2 [^ (r) - ^ (-r) + & (-r) - ^ (r)] = 0. (46) 

By the boundary conditions on the potentials 

1>j ( r ) ~ 1>i (- r ) < oo. (47) 

That is, it is finite. Therefore ipj (r) — ipj (— r) + ipi (— r) — ^ (r) is also finite. 
According to the theory of harmonic analysis, any finite function satisfying the 
Laplace equation is a constant. Therefore, we may set 

ipj (r) - ipj (-r) + i[>i (-r) - ipi (r) = 
15 



or 

^ (r) - ^ (-r) = fa (r) - ^ (-r) . (48) 

That is, the odd parts of the potentials are equal independently of i and j. This 
means we may set 

rl>j (r) - ^ (-r) = ft (r) - ^ (-r) = 2y (r) , (49) 

where Y (r) is the odd part of the potential. 

The symmetry conditions for the even part of the potential is obtained as 
follows. From Eq. (ITS1) taken to O (e 2 ) 

- V • w« [ez,n 2 V^ (r) + fc s TV/ l2 (r)] = 0. 

On integrating it over $7 and observing the boundary conditions discussed ear- 
lier, we find 

= V [-uj^n 2 [fa (r) + ^ (-r)] - 2k B Tf u (r)] (50) 

which means 

-uJiezin 2 [ipi (r) + fa (-r)] - 2uj l k B Tf il (r) = c, 
but ipi (oo) = and /« (oo) = n 2 . Therefore c = n 2 and we find 

fa (r) = /« (-r) = n 2 - |g^ [^ (r) + fr (-r)] (51) 

and similarly 



,2 



/tf W = fn (">") = n' - ^ [^ (r) + ^ (-r)] . (52) 
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Upon use of the Poisson equations (|4"4")) an d P5|) , it follows 



V 2 [^ (r) + 4> 3 (-r) + ^ (r) + ^ (-r)] 



'" ' ' : ' ^ ; H [^ (r) + ^ (-r) + fr (r) + fr (-r)] , (53) 



Dk B T 



(54) 



which yields 

( v ' v_ y) ^ (r) + ^ (_r) + ^ (r) + ^ ; ( ~ r)] = °- 

Note that z = |zi| = \zj\. By the boundary conditions 

[tpj (r) + V'j (-r) + tj>i (r) + -0. 4 (-r)] = finite 

or 

i\)j (r) + i\)j (-r) + Vi (r) + Vi (-r) = 0. (55) 

Therefore there exists a function T (r) independently of i and j such that 

^ (r) + ^ (-r) - - [^ (r) + ^ (-r)] = 2r (r) . (56) 

Here T (r) is the even part of the potential. Thus we have now identified the 
even and odd part of the potential. Therefore we finally obtain 

^(r) = -Vi(-r)=r(r)+y(r), (57) 

^ (-r) = -ft (r) = T (r) - Y (r) . (58) 

These are symmetry conditions for ionic potentials. 
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2. Diagonal Distribution Functions 

With conditions (1571) and ([55]) we deduce the diagonal parts of the distribu- 
tion function matrix 

/« (r) - fa (-r) - n 2 + §^ [^ (r) + fr (-r)] = n 2 + f^r (r) , (59) 
fii (r) = /,, (-r) = n 2 - Z £^ fa (r) + ^ (-r)] = n 2 - |^r (r) . (60) 

3. Cross Distribution Functions 



It is expected from Eq. (129)) that cross distribution functions also consist of 
even and odd parts. We will write them as 

fa (Tr) = fji (±r) = n 2 + G(r)±U (r) . (61) 

Before proceeding further we would like to note that in Harned and OwerP^J n 2 
is inserted in the even part whereas in Wilson's thesis the factor n 2 is absent. 
If fij ~ n 2 as r — >• oo, then the boundary value of G(r) would be G(r) — > n 2 . If 
fij is a perturbation to the equilibrium distribution function as it is here, there 
is no need for addition of n 2 . We have put n 2 in Eq. (|6ip assuming it is for the 
entire distribution function. 

D. Differential Equations for Even and Odd Parts 

Inserting Eq. (|61[) into Eq. (|29p. separating the symmetric and antisym- 
metric parts, and making use of the Poisson equations, we finally obtain four 
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coupled second-order differential equations for G, U, T, and Y: 



72 1 2\-/x !j„2> w _ ,.,9U(r) 



where 



V-- - -^ j G(r) - -rfn^T (r) = -/*'-^, (62) 

V--^)^)-^, (63) 

V 2 r(r)-^r(r) = ^ 7 G(r), (64) 

V2y(r) = W CI(r) ' (65) 



" = teT' (66) 



with 



for binary electrolytes. Solutions of the equations given above would yield the 
potentials and distribution functions with which to calculate average quantities 
for the binary electrolyte. 

IV. Formal Solutions 

Owing to the fact that the external field is applied in the positive x direction 
the system is axially symmetric. Therefore we will adopt the cylindrical coordi- 
nate system with the x direction as the axis of the cylinder. The coordinates are 
denoted (x,p,9), where 9 is the azimuthal angle and p is the radial coordinate. 
In this coordinates the Laplacian is given by 

V 2 -— -— — — — 
dx 2 p dp dp p 2 d9 2 

Because of the axial symmetry, the solutions are independent of angle 6. The 
angular derivative in V 2 therefore can be ignored. The functions G, U, T, and 
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Y then depend on x and p only. We define one-dimensional Fourier transforms 
of these functions as follows: 

G(x, P ,0) \ 2 f°° ( g(a,p) \ .... 

V% P ^) = ir das ^ {ax) (u(a,p)\ 
Y (x,p,0) J tt J V 2/(«)P) / 

where a is the Fourier transform variable. It should be noted that odd functions 

(U, Y) are sine transforms whereas the even functions are cosine transforms. 

Taking Fourier transforms of Eqs. (|62j) -(f65 |) according to Eqs. (|69| and (f70|) . 

we obtain differential equations for g (ap) , 7 (ap) , u {cap) , and y (ap) : 

-— p— -a 2 - — ) u(a,p) = p.'ag(a,p), (71) 

p dp dp 2 J 

1 d d k 2 \ 1 

-^-P^--a 2 - y J 9(a,p)- ^v'n 2 K 2 -f{a,p) = -p/au(a,p), (72) 

-^-PTT "« 2 ^ -tt) 7(«>P) = ^-^-7ff(a,P), (73) 
p dp op 2 / 2n^?7' 

-l-P^ - -" 2 ) y(«>P) = T~2~7"( a '^)- ( 74 ) 
p dp dp J 2n z r]' 

The boundary conditions on potentials are deduced from the boundary condi- 
tions already established: 

lim P 5^1 = 0. (76) 

p-s-0 C»p 

This set of equations can be solved more concisely if an equivalent eigenvalue 
problem is solved for the set. However, to remain as close as possible to the 
original approach, we follow his approach. 

Eliminating u(a, p) and 7 (a, p) between Eqs. (|TTj) . (jT2"j) . and (1751) we obtain 
a fourth-order differential equation 

-^-P^--c 2 J g(ap)=[-K i -p' 2 a 2 \g(a,p), (77) 

c 2 = a 2 + \k 2 . (78) 
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The solution of the homogeneous differential equation of Eq. ([77]) 

(&£-*)*-* <79) 

which is finite at p = oo, is found to be the zeroth-order Bessel function^ of 
second kind K (cp) : 

4>{p)=(j> K {cp), (80) 

Therefore the solution of Eq. (|77|) . g(a,p), must be a linear combination of 
zeroth-order Bessel functions of second kind. The coefficients arc determined as 
follows: Set 

g(a,p)=w(a)K (\p), (81) 

where ui (a) and A are determined such that ui (a) Kq (Xp) is a solution of Eq. 
([77]). On substituting Eq. (JHJ) into Eqs. (J77J) we find 



( \1 2\ 2 ( 1 4 /2 2 \ 

(A -c ) = I -« - /u a I 



and the solution for A 2 is 



A 2 = c 2 + ]-k 2 R = a 2 + |k 2 + ]-n 2 R, (82) 

A 2 = c 2 - i« 2 i? = a 2 + ]-n 2 - \k 2 R, (83) 

where 

The existence of two values of A implies that the solution for g(a, p) is the linear 
combination 

g(a, p) = Wi (a) K (Xip) + w 2 (a) iT (A 2 p) . (85) 

Similarly, the functions u(a,p), 7 (a, p), and y (a, p) can be determined as the 
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linear combinations 



7 [a, P) = Xi (<*) I<o (Aip) + X2 (") ^o (A 2 p) + X3 (a) ^o (A 3 p) , (86) 

u (a, p) = ipi (a) X (Aip) + tp 2 (a) K (X 2 p) + ~ip 3 (a) K (A 3 p) , (87) 
V (a, z 5 ) = £i (a) #b (Aip) + £,2 (a) K (X 2 p) + £ 3 (a) ^o (A 3 p) 

+ £ 4 (a) # (A 4 p) , (88) 



where 



\ Z Z Z , Z 

X^ = c = a + -k 



A4 = a. 



(89) 
(90) 



The coefficients uj\ (a), \i ( a )> $1 (<*)> £1 (a), e ^ c - are determined on substi- 
tution of Eqs. (JHSJ)— dHHl) into Eqs. ^ -flU]). Since i^ (Aip) and K (A 2/ o), 
■Ko (A 3 p), and -Ko (A 4 p) are independent, the aforementioned substitution yields 
the set of equations for the coefficients: 

1 



\2 „2 



—k Rtpi {a. 
-K 2 Rip 2 (a 
— J wi (a) + -Tl'n 2 K 2 xi (a 



A9 — a z "— ) W2 (a) + T.il'n z K' i X2 (a 



I/22 
—?i r • - - 

2 
-r]'n 2 K 2 X3 (a 



\2 „2 



— /i'awi (a) = 0, 
+ /j,'au!2 (a) = 0, 
+ /i'a-01 (a) = 0, 

+ ^'cm/>2 (a) = 0, 

- p,'aip 3 (a) = 0, 
0, 



K 2 \ K 2 

2 / 2n^?7' 

k 2 \ k 2 

ir X2 (a) ~ 77^-7^2 (a) 
2 / zn^r]' 



\z n z 
\ 2 -a 



0, 



(A 2 - a 2 ) 6 (a) - ^^r (a) = 
(X 2 2 -a 2 )t2(a)~^^(a) = 
(A! - a 2 ) 6 (a) - ^> (a) = 



(91) 
(92) 
(93) 
(94) 
(95) 
(96) 

(97) 

(98) 

(99) 

(100) 
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From the boundary conditions follow the equations 

3 



^2xi (a) = 
1=1 


= D 


!>(«) = 


= 0, 



(101) 
(102) 



i=i 



while since the flow is sourceless and the field is divergenceless it follows 

0. 



lim p 

p->0 



2 ,9j{a,p) dg{a,p) 
n r\ 



dp dp 

llm n P A =0 ' 

p^o dp 



which provide additional equations 



n V (Xi +X2 + X3) + wi + w 2 = 0, 

ipi + i>2 + i>3 = o. 



Therefore on combining Eqs. (J10ip - (|104|) we obtain 



Ul\ + CJ2 = 



Xl + X2 + X3 



2 / 

en 7] 



D 



(103) 
(104) 



(105) 
(106) 



There are 13 equations for 12 variables. These equations are first reduced to 
equations for Xi only. Then all other variables are determined from \i- We 
obtain 



Xi(«) 
Xa(a) 
Xs(«) 



e(l + fl) 
2L>i? 2 ' 

e(l--R) 
2DR 2 ' 



1-i? 2 



Z?i? 2: 



(107) 
(108) 
(109) 
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which give rise to the solutions 



uj\ (a) 


u> 2 (a) 


4>i («) 


ip 2 (a) 


ip3 (a) 


6 (a) 


6(«) 


6(a) 



2 ,e(l + i?) 
' 2DR 
2 ,e(l-R) 



-n r\ 



2DR ' 



9 , Vl — i? 2 e (1 + i?) 

n V — K~ 



DR 2 



2 .Vl-R 2 e(l-R) 
1 2 DR 2 



2 , 2a// e 2 / A m e 

Vl-1? 2 " e 



2 Di? 2 ' 

VI - R 2 e 

2 Di? 2: 
e 



Vl-i? 2 



£i? 2 ' 



6=0. 



(110) 

(111) 

(112) 
(113) 
(114) 
(115) 

(116) 
(117) 
(118) 



Finally, the solutions for g[a,p), etc. are obtained: 
'(a»p) 



w (a, p) 



n 2 ej \{l + R) 

D 
n erf n 



Dh 



2R 

a(l + R) 



Ko (Aip) - ^^-K (A 2 p) 



27? 



(119) 



/i>:i ^o (A lP ) + a{1 R2 R) K (X 2 p) - gtf (A 8 p) 



7(a.P) = -j 



! ' R) ^o (Xip) + KJ^K, (\ 2P ) l —^K, (X 3 p) 



2R 2 



2R 2 



R 2 



y(oi,p) 



ep 



q/ q/ 2a 

— K [Xip) + -j^Ko (X 2 p) - -jp K v (Xzp) 



(120) 

(121) 
(122) 



Dk 2 

with which the Fourier transforms of the final solutions for the distribution 
functions and potentials can be constructed in the manner of Eqs. (|69[) and 
(|70[) . For this purpose it is convenient to define the following abbreviations for 
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integrals: 


Gi(p) 


G 2 ( P ) 


Ui( P ) 


u 2 ( P ) 


u 3 (p) 


ri(p) 


t 2 (p) 


t 3 (p) 


n( 


Y 2 ( 


Y 3 ( 



/I _i_ D\ 

da — — — cos (ax) K (Aip) , 

(1 — Ft) 
da — — - — cos (ax) K (X 2 p) , 
IK 



2 r a(l + R) 

da —5 sin (ax) K (MP) , 

R z 

da —5 sin (ax) K (X 2 p) , 

R z 

da— sin (ax) K (X 3 p) , 



2 f°° (\ 4- R) 

da — ^-^ — cos (ax) Kq (Xip) . 



7T 







2R 2 



2 roo /1 r>\ 

~ / da on9 C0S ( aa; ) ^0 (A2/9) , 
7T 7 2R Z 

1 r>2 

da — — ; — cos (ax) Kq (X 3 p) , 

7T./0 « 



2 f°° a 

p) = — / <ia — ^ sin (ax) Kq (Xip) . 
7T Jo R 

2 f°° a 
p) = — da—r sin (ax) Kq (X 2 p) . 

71" Jo R 



P) = - 



2a 

da—; sin (ax) Kq (X 3 p) . 
R z 



(123) 
(124) 

(125) 
(126) 
(127) 

(128) 
(129) 
(130) 

(131) 
(132) 
(133) 



The formal solutions for the distribution functions and potentials are given by 
the forms: 



2^2 






fjj W = n 2 



2^,2 



Zje n 
Dk B T 



[Ti ( P ) + r 2 ( P ) + r 3 ( P )} , 
[Ti (p) + r 2 ( P ) + r 3 ( P )} , 



ho (Tr) - /ji (±r) 



[Gi (p) + G 2 (p)] 



2 / 
n evj 

D 

±^f[U 1 (p) + U 2 (p) + U 3 (p)} 



(134) 
(135) 

(136) 
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and 



i>i (r) = -tpi (-r) 



J p! ( P ) + r 2 ( P ) + r 3 ( P )] + ^ [Ki ( P ) + r 2 ( P ) + r 3 ( P )] , (137) 



V-i (~r) = -ipi (r) 



-^(pj + raOoj + rafp)] 



e//, 



[Vi (p) + Y 2 (p) + Y 3 (p)} . (138) 



These formal solutions given in terms of quadratures of the zeroth-order Bessel 
functions of second kind can be used to calculate various local mean quantities 
associated with transport in binary electrolyte solutions. Since the Bessel func- 
tions are not a trivial function to compute the quadratures with, it is useful to 
review briefly their mathematical properties we can make use of to evaluate the 
quadratures in question. 

The Bessel function K v (z) obeys the differential equatio n 1 23 ! 24 ! 



1 d d 

z dz dz 



,2\ 1 



1 



K v (2) = 0. 



(139) 



It is a regular function in complex z plane cut from z = to z = - oo. One of 
its integral representation is given by 



/■oo 

K v (z)= dte- zcosht cosh(vt) (|argz|<- 

Jo V 2 



(140) 



We will make use of this integral representation. We also note one of the recur- 
rence relations 

K' v (z) = -K y+1 (z) + -K v (z) . (141) 

z 

Here the prime denotes differentiation with respect to z. This recurrence relation 
implies 

K' (z) = -K, (z) . (142) 

In particular, K (z) may be represented by the ascending series 



K (z) 



ln(-)+ 7 



^oW + Ef 1 



i=l 



1 



1\ (I)' 



(ny 



(143) 



2G 



where Iq{z) is the other independent solution of Eq. (|139[) , which is given by 
the ascending series 



oo I z 



, 27! 



m^e|f- ( 144 ) 

This function is an entire function if the cut z plane. The asymptotic behaviors 
of Kq (z) and Iq(z) are follows: as \z\ — > oo 

K (z) = J^e-*[l + O(z- 1 )], (145) 



h(z) = ^/[l + 0{z- 1 )} (146) 

for |argz| < f . It is useful to note that the presence of In (|) in Ko (z) clearly 
indicates that Kq (z) is not only discontinuous across the branch cut, but also 
diverges logarithmically as z — >• 0. If Kq (z) is analytically continued across the 
cut, we obtain, for example, 

K (e iw z) = K (z) - ml (z) . (147) 

This relation is easily seen to be true from Eq. f)143|) . 

Wilson makes use of only Eqs. (|142[) and (|143|) to evaluate integrals for the 
velocity and ionic field, which involve the Bessel functions. This analytically 
continued form as well as the properties presented above can be made use of for 
exact evaluation of the quadratures presented in Eqs. (|123|) (jl33[) . as will be 
shown elsewhereP^S 

V. Electrophoresis 

Onsager in his 1927 papei^ proposed that there exist two effects in elec- 
trolytic conduction, one the electrophoretic effect and the other the relaxation 
time effect. In Ref.Pthe electrophoretic effect was treated by making use of the 
Stokes law^ and the relaxation time effect by estimating relaxation time by a 
heuristic argumentPS first used by J. C. Maxwell. In the Onsager- Wilson theory 
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the electrophoresis effect is calculated by solving the Navier-Stokes equations^ 
and the time of relaxation effect by calculating the ionic held using the solutions 
presented in the previous section. In this section, we discuss the electrophoresis 
effect. 

The Navier-Stokes equations^ necessary for the calculation of the effect 
may be cast into the forms 

770V x V x v=-Vp + F, (148) 

V-v = 0, (149) 

where F is the local external force, p is the pressure, and 770 is the shear viscosity 
of the fluid. In solving these equations Wilson makes the following assumption: 
the forces exerted by the field due to an ion and the ion atmosphere upon the 
charges in the ion atmosphere may be neglected, because they are proportional 
to the square of the charges of the central ion. This means the external forces 
include forces linearly proportional to the external force X and the ionic poten- 
tials. Therefore 

F=pX 

and since the charge density j is described by the Poisson equation 

n v 

F x = ——V% (r) (150) 

and the Navier-Stokes equation may be put in the form 

770V x V x v = -Vjj + F (151) 

for a binary electrolyte. 

To solve Eq. (|151|) we look for a vector a satisfying 

7/0V 2 (V 2 a) = F, (152) 



2S 



where F is the vector form of Eq. f|150[) . Eliminating F between Eq. (|148[) and 
Eq. CHU) , we obtain 

ry V x V x v + Vp - r/ V 2 (V 2 a) = 0. (153) 

Noting the identity in vector algebra 26 

V x V x (V 2 a) = V [V 2 (V • a)] - V 2 (V 2 a) , 

it is possible to express Eq. (J153I) as 

?7oV x V (v+V 2 a) +V [p - ??oV 2 (V • a)] = 0, (154) 

but 
V x V x (V 2 a) = V x V x [V 2 a-V (V • a)] = -V x V x V x V x a (155) 

because 

V x V x V (V ■ a) = 0. 

Inserting Eq. (|155|) into Eq. (|154|) . we obtain 

r/ V x V x (v-V x V x a) + V [p - rj V 2 (V • a)] = 0. (156) 

Since two terms in Eq. (|154[) are independent, it is found 

v = VxVxa, (157) 

P = Po + %V 2 (V • a) . (158) 

In Eq. (|157p for v, the sign error in Wilson's thesis is corrected, and r/o is also 
added in Eq. (|158[) . In any case, the solution of Eq. (I152|) provides the velocity 
field and the local pressure arising from the local potentials. It, however, should 
be noted that the solution of the equation 

V [p - 7? V 2 (V ■ a)] = 

29 



is more generally 

p-r] V 2 (V -a) = p , 

where po must be uniform in space in the hydrodynamic scale for the Navier- 
Stokes equation, so that 

Vp = 0. 

Therefore, po may be regarded as a uniform equilibrium pressure in the absence 
of the external field. Thus, to be correct, po is added to Eq. (|158p . The po is 
missing in Wilson's expression (1158[) for p. Here po is the equilibrium pressure 
of the solution that may be given, for example, by the virial form of pressure 20 : 

po = nk B T -^Yl n ^ r *"*^r-9% M ' ( 159 ) 

where Uij is the intermolecular potential of pair (i,j) and g . (r) is the equilib- 
rium pair correlation function. The index i runs over the species in the fluid, 
including ions. Determination of g°- (r) should be made by following the modern 
theory of equilibrium Coulomb (ionic) fluidsp2l 

It is now necessary to find the vector a by solving Eq. (|152[) . Since the local 
force F x is a linear combination of the solutions for potentials, it may be written 
in the form 

3 3 

F x (z, p) = m J2 CiK (Xi P ) +VoY, S ^« ( A ^) ' ( 16 °) 
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where G; and Si are integral operators denned in a space as below: 
„ eX [°° ^ (l + i?)(A 2 -q 2 ) 

eX f°° , s (1 - R) (A2 - a 2 ) , . 

C 2 = --^— / da cos (ax) -^ S (161) 

^ Vo Jo 2R 

^ r , , A 1 - r2 ) ( a i - « 2 ) 

C 3 = , / dacos(ax)^ ^ S 

eXy' [™ . , a(A 2 -a 2 ) 
5,1 = ~o~~ 2 2 / "asm (ax) 



27r 2 77 K 2 J Q y ' R 2 

eX// Z" 00 2a (Af — 



eXyi [°° . a (A 2 - a 2 ) 



S 3 = —^ j / da sin (ax) 

Z7T ryot Jo 



R 2 

With i 7 ^ so expressed, the equation for the x component of vector a, Eq. (|152[) . 
reads 

3 3 

V 2 (V 2 a x ) = Y, C t K (V) + J2 S * K o (Aip) • (163) 

i=l i=l 

This is solved by 

v a o a =y;(c i +s i )§^+>i* ) (164) 

where A* is the solution of the homogeneous equation 

V 2 (V 2 A*) = 0. (165) 

It must be also a linear combination of Kg (\ip), but since it must satisfy the 
boundary conditions we choose it in the form 

A-— E(C* + Si)#^. (166) 

»=i * 



Therefore we have 

V a x = ) (Cj + Sj) -2 5 . (167) 

*-^ \f — or 

i—l l 

It should be remarked that this solution is opposite in sign to Wilson's. The 
solution of this equation must be also a linear combination of the Bessel func- 
tions Kg (Xip) and Ko(ap). Wilson does not make the solution process for this 
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equation explicit. Here we make it transparent. The solution is sought in the 

form 

3 1 
a x = V {d + Si) - 2 - {h [Ka (X iP ) - K (ap)] + b 2 [K (ftp) - K a (ap)}} , 

% — \ l 

(168) 
where j3 is a parameter to be determined such that the right hand side is a 
solution of Eq. (|167[) . On inserting this into Eq. (|167p we find 

[h (A? - a 2 ) - 1] K (Xip) + b 2 (ft 2 - a 2 ) K (fiip) + K (ap) = 0. 

Now, b\ and bi are chosen such that 

bl = ' (169) 



and 



that is, 



lim [6a (ft 2 - a 2 ) K ((3 lP ) + K (ap)] = 0, 



lim b 2 (ft 2 - a 2 ) Ko{(3 lP ) = -K (ap). (170) 

Pi — y& 



Then we find 



Now, since 



1 
{P 2 ~a 2 Y 



b2 = -7^—T V (171) 



r 1 ta \ r ( \\ l r [^o CM - K (ap)} 

llm 7132 2\ i K ° (PiP) - K ® ( a P)l = T~ llm n5~, 1 

fo^ot (ft 2 - a 2 ) r/ r/J 2a p^a (ft + a) 

= £*S(«p) 

= -£^i(«p) (172) 

in which we have used recurrence relation (|142|) for the last equality and the 
prime denotes the derivative with respect to the argument ap, we finally obtain 
the solution of Eq. (|167|) in the form 

-X^m c ^ / [ X ° i^P} ~ ^° ( a ^)] ! P Kl ( a P) 1 ri7Q^ 

^-^(Q + SOj ^—^ + 2a(Af _ a2) j- ( 173 ) 
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With this result for a x it is possible to calculate the velocity by using the 
expressions 

clival ]T [_Ci + St] \ - [K (X z p) - K (ap)} + P K^ap) 

(174) 

V, diva= - J2 [Ci + Si] I ^ 2 LMV) - ifo(ap)] + / ^M 

(175) 

A. Axial Velocity 

Inserting the formula for a x into Eq. (|157[) and Eq. (|158[) and making use 
of Eqs. ()174j) and (|175|) we readily obtain the velocity field. The axial velocity 
is finally obtained as below: 

v x (x, p, 0) = — -= — / da cos (ax) x 
27H770 Jo 

" /,,)A L[K (A 2 p)-if (^)] 



2R 2 (Al - 

1 ' " 7/ "' A ' [^o(A 3 p) - K (ap)] + ^lfi(ap) 



i?2(A|-a 2 ) 1 " uv ' ,<,Fy " UV "" J ' 2 

el/j' f 00 f aA? 

da sin (ax) <^ 2 — [K (X 1 p) - K Q (ap)\ 



2TT 2 r] K 2 J y '{R 2 (X 2 -a 2 ) 

+ R 2 (Xl-a 2 ) [Ko{X ^- K ^ ap)] 
2aX 2 



R 2 (X 2 -a 2 ) 



-, -^ ,K (X 3 p) - K (ap)} <► . (176) 
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in which the integral operators Cj and Sj defined in Eqs. f|161[) and (|162j) are 
made explicit. Further simplifying it, we obtain 

V~ /"OO 

V;e (x, p, 0) = — — — / da cos (ax) x 
2tH?7o Jo 

[A^o(Aip) + X 2 2 K (X 2 p) - 2A 2 (1 - R 2 ) K (X 3 p)] 



K 2 R 2 



-^o(ap) + -y 



' dasin(ax)^- I [A^o(Aip) + A2x (A 2 p)-2A^ (A 3 p)] 



2ir 2 r]oK, 2 J k 2 R 2 

(177) 

Since the quantity of interest for electrophoresis is the axial velocity at the 
center ion (the coordinate origin) Wilson at this point consider the special case 
of v-c (x, p, 0) at x = and p — before fully evaluating the integrals in Eq. 
(|176|) . That is, if x = and p = are taken, the axial velocity component at 
the center ion is given by 

eX f°° , f A?(l + i?) , /Ai 

A|(l-fl) , /A 2 \ A|(l-i? 2 ) /A 3 \ l\ 
" 2/P (A 2 - a 2 ) ln UJ + fl 2 (AI - a 2 ) * UJ ~ 2 J ' (1?8) 

He then evaluates these integrals by using the method of contour in complex a 
plane. For the purpose it should be noted that 



A? (14 


R) 

-a 2 ) 
R) 

-a 2 ) 


— 


K 2 


1 + t z + Vl - 2x 2 t 2 


2R 2 (Xj 


2 

K 2 


1 - 2x 2 t 2 


Ag(l- 


l + t z -Vl- 2xH 2 


2R 2 (X 2 


2 


1 - 2.T 2 * 2 


A|(l- 


R 2 ) 


= 


x 2 


2t 2 (1 + 1 2 ) 


R 2 (X 2 


-a 2 ) 


\-2xH 2 ' 




x = 


K 




K 



where 



These expressions mean that there are simple poles at t = ±1/ (y2x) in the 
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integrands. The arguments of the logarithmic functions have branch points 



t = ±i v / 2(l + x 2 ) for Ai 
t = 0, — oo for A2 
t = ±i for A3. 

Furthermore, the integrands of the first three integrals vanish to zero as \t\ — > 00. 
Therefore the first three integrals can be evaluated by applying the method of 
contour integration. Wilson asserts that the integral of — | for the last term can 
be replaced by a contour integral C used for the other integrals and the contour 
integral vanishes because the integrand is constant and hence 

Thus the last integral does not contribute to the electrophoretic effect. On this 
ground, he obtains from Eq. (|178[) the formula for the electrophoretic effect as 
represented by the formula 

„ x (0,0,0) = -^^-/(x) (179) 

where 

fix) = H =- (2a; 2 sinlT 1 x + V2x - x\/l + x 2 

- (1 + 2x 2 ) tan" 1 (V2x) + (l + 2x 2 ) tan" 1 X 1 . (180) 

\ / v l + x 1 ) 

In effect, this result is correct, provided the last integral does not contribute 

to the electrophoretic effect or negligible. However, as it stands, Eq. (|1T8[) 

definitely contains a divergent integral. 

In fact, in the case of x — the integral in Eq. (|176[) 



I da—Ki{ap) 
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can be exactly evaluated for all values of p > 0. By using K' (z) = —K\(z) and 
the integral representation for Ko(ap) it is easy to show that 



^c 



daaK^ap) = — , (181) 

o Z P 

which indicates the integral is singular at p = 0. Therefore p — should not be 

taken within the integral. In fact, if the ions are assumed to be hard spheres as 

is usually assumed in the treatment of electrolytes — in fact, this is what is done 

in the Debye-Huckel theory underlying the Onsager theory of conductivity — it 

may be appropriate to take the mean hard sphere radius, say, a for the minimum 

value of p; p = a or even p = n^ 1 , the radius of the ion atmosphere. Then the 

axial velocity emerges finite for a < p < oo. But then, it is necessary to evaluate, 

for example, the integrals for v x (0, a, 0) with the form 

v, (0, a, 0) = -^- f°° da I -11±^4- [KoiXra) K (aa)} 
2?r Vo Jo l 2R (Ai - or) 

(l-^)Al 



2R 2 {XI - a 



[K {X 2 a) - K {aa)} 



i 1 -* 2 )^ [KoiX^-KoiaaA-^L. (182) 



R 2 (X 2 3 -a 2 ) L uv ° ' uv /J J 87r7?o<r 

The integrals in this expression, however, cannot be evaluated by the method 
used by Wilson. It would be straightforward to evaluate the integrals numer- 
ically on computer. There then would be an extra contribution to the elec- 
trophoretic effect arising from the last term on the right in Eq. (|182[) . and as a 
matter of fact, the terms associated with Ko(aa) in Eq. (I182[) all together also 

contribute 

eXn 



Sirrjo (np) 



to the velocity. It also diverges as p — > 0. If these are evaluated at p = a 
or p = k^ 1 , they are not negligible at all and hence would modify Wilson's 
electrophoretic effect significantly. 
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It is therefore possible to conclude that Wilson's result for the electrophoretic 
effect represents only the nondivergent part with the divergent contribution 
arbitrarily discarded. This divergence difficulty might have been the underlying 
reason for not publishing the one-dimensional Fourier transform approach taken 
in Wilson's treatment in favor of a full three-dimensional Fourier transform 
approach that had to wait twenty years until the Onsager-Kim theory^ see 
page 146, section 8, Chapter 4 of Ref.pS 

VI. Ionic Field 

To calculate diffusion of ions it is necessary to know the ionic field of ion j 
located at x = — the center of ion atmosphere. More precisely, we need the 
contribution to the ionic field arising from the ion atmosphere. It is given by 



AX (r) = - W,- (r) . 



(183) 



Since ipj ( r ) is already known it is trivial to compute. Inserting the formula for 
4>j (r) obtained earlier we find it in the form 



AX(r)= V / da 

kD J 



cos (ax) 



{1 + r) k (x ip)+ { ±-Vk (x 2P ) 



2R 2 

I -R 2 



2R 2 



R 2 



K (\ 3 p) 



TV 



2e/i' 



r o. ol 

da sin (ax) [~^Ko (Aip) + ~j^ K o (Mp) (184) 

--j^K (\ 3 p) 
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and hence 



2e f°° 

AX (r) = e x / daa sin (ax) 

kD J 

. (1--R) 

2B? 

T e x — -pr^ / daa cos (ax) 
■kDk 1 J 

2e ^ r , , x 

aa cos (ax) x 



2i? 2 



#o (Aip) 



K (X 2 p) - l -^~K (X 3P ) 



-^i^o (Aip) + ^2-Kb (A 2 p) - -^i^o (A 3 p) 



^7T^ 2 i 



(1 + f 2 )Al ^(A lP ) 

2i?2 IV 1AM 



{1-R)\ 2 „ ,. , , (1 - fl 2 ) A 3 ^ 

-ifi (A 2 p) + -^ — Ki {X 3 p) 



2B? 



2e ^ r% . / x 

T e p _ n „2 / eta Sin (ax) 



'ttL>k 2 



R 2 



(185) 



2aAj 

~R2" 



Kx (X 3P ) 



Therefore if we set x = then 



AX(0,p,0) = e x AX x (0,p,0) 



Te x - 



2eu' r x . a 2 



da-j [K (Aip) + X (A 2 p) - 2X (A 3 p)] • (186) 



Wilson calculates the force on ion j in the axial direction by further setting 
p = in this equation, which then is given by the expression 



AX(0,0,0)=± e;c 



ep' 



a 2 , (X\X\ 



„ ^ ln V A 



He then evaluates it by integration by parts in the form 
ep! k 



AA(0,0,fl)=T 

g(x) = - 



2D 
1 

2^3 



■ff(^), 



(187) 



(188) 



■\/l~ 



x + tan 



ViT^ 5 



V2x - tan -1 [v^Xj 
(189) 



The original formula in Wilson's thesis is different from this, but this is the 
formula in the summary of his work given in the monograph^ of Harned and 
Owen. We may verify it by a different integration method as follows. The 
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integral in Eq. (|187l) may be written as 

v2 /\2\2' 



Il = 



, or , / AfAa 
dt- 



72/ 2^,/n 



i 2 



[in A? + In A;;- 2 In A3] 



Using the properties of Xi given earlier for the axial velocity [see Eqs. 

and Eqs. (155)) ] by applying the method of contour integrations^ similar to those 

for axial velocity integrals, we obtain 



Il = 



AV2e 

1TK 3 

4\/2£ 2 

7TK 3 
~4^" 



dt-z— 



V 2 (!+« 2 ) 



dj/ 



a 2 + y 2 



(-7U) + 2 



- 2 / dy- 
'o 



f/f- 



f 2 



(-7T*) 



(t = complex) 



a 2 + y 2 
tVO- + H + V^ + ^ tan" 1 (2CV1 + C 2 ) - tan" 1 (V2£ 



(190) 



Here we note that 



~ tan- 1 fa^/TTe) = tan- 1 -= 



by the addition theorem of tangents. On inserting this integral into Eq. (|187[) , 
we obtain g (£) as in Eq. (|189[) . The result for g (£) in Eq. f|190[) agrees with 
Wilson's. 

This completes the critical review of the calculations of the axial velocity and 
ionic field and Wilson's evaluation of /(£) and g(£) on the basis of Onsager's 
theory of conductivity. We have pointed out that taking p = within the 
integrals before they are evaluated can cause a difficulty because of the lack 
of uniform convergence of some of the integrals. If p ^ then an alternative 
method of evaluating integrals must be used. An alternative method will be 
discussed elsewhere^! that yields exact analytical results for the integrals for 
all values of p. This method recovers Wilson's results as the leading order 
contributions and the correction terms are at least O (n 2 p 2 ). 



39 



In any case, the axial velocity formula (I176P given here does not appear in 
the existing literature. Our analysis of the axial velocity will start from Eq. 
(|176[) in the sequef^D to this tutorial review article. 

VII. Conductance 

The total field strength E is given by 

E = X + AX :J (0), (191) 

where AX,(0) is the field arising from the ion atmosphere interacting with X. 
This term is called the relaxation effect, but in effect it is a local dressed field 
in the modern terminology. The velocity Vj of ion j in the effective (dressed) 
field is then given by 



„x(l + ^)) + v,(0,0,0). 



(192) 



Here v^ (0, 0, 0) is Wilson's axial velocity. It is useful to remark that the diffusion 
current (velocity) may be written relative to v x (0, 0, 0) as^ 

Ji=Vj-v a (0,0,0). (193) 

Then it is seen that the first term on the right of Eq. (|193[) is the diffusion 
flux in the absence of a density gradient and that the so-called electrophoresis 
effect is nothing but the reference velocity suitably chosen in linear irreversible 
thermodynamic theory, or kinetic theory, of transport processes. In fact, the 
calculation of w x (0, 0, 0) clearly supports this interpretation. Thus inserting the 
results for AXj(0) and v x (0, 0, 0) we obtain 

Ce^Uju'K led ft \ 

W-^rfO-jjjTjU/tt)). (194) 

The mobility in electrostatic units is 



X 

40 



(195) 



and in practical units 



1 / e^ojj/i k 



W« m 



At infinite dilution k — > 0, and the limiting mobility is given by 



300 



Since the limiting conductance is 



A° = 964931u°, 



we have 



A - A - ^ 
J ~ 3 Dk B T" 3 



A?ff(0 



6 x 300^7^0 
With the definition of equivalent conductance 

A = A + + A_ = Ai +A 2 , 



mi 



(196) 



(197) 



(198) 



(199) 



the equivalent conductance for the electrolyte is given by 
A(0 = A°- 



Dk B T 9{t;> 6 x 300^^0 /lU 



In experiment, the relative equivalent conductance 



AA(0=A(0-A(0) 



(200) 



(201) 



is measured and reported. We use the ratio of the relative equivalent conduc- 
tance 



AA(0 Z>4fAV0) 



1- 



3(0 
9(0) 



(|ei|+|e 2 |)K r/Q\ 



1 - 



M) 

/(o) 



(202) 



A(0) A(0) 

This expression can be used to compute conductance and compare it with ex- 
perimental data. 
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VIII. Concluding Remarks 

Charged particles in condensed phase in an external electric field are com- 
mon occurrence in physics, chemistry, and biology. It, especially, is important to 
comprehend theories of transport properties of such systems and have practical 
theoretical tools for them from statistical mechanics standpoints. Among the 
existing theories, the theoretical development initiated a long time ago by On- 
sager and his collaborator j 9 l 10 l 1:L l 12 l 18 l provided useful theoretical models that 
after complex mathematical treatments have yielded some useful well known 
analytical results, but his theories have not received critical reviews by other 
workers that might have stimulated studies for possible further developments in 
the field and other related fields, which are many currently. His theories should 
be as relevant today as they were useful when they were developed. Given the 
power of Onsager's analytical theories on the subject of electrochemistry, the 
apparent scarcity of critical follow-up studies of the theories^SES is quite curi- 
ous. The present author sees that there seem to be considerable insights and 
lessons one can draw from them only if we critically examine the theories and 
benefit from them for the current topics of interest. Motivated by this thought 
and their possible utilities of the theories for the currently investigated phenom- 
ena in charged condensed matters subjected to external electromagnetic fields, 
the present author has been studying his theories. In this tutorial review ar- 
ticle one of his early work described in Wilson's PhD thesis has been studied 
critically and in detail. It is found that it has a difficulty in one respect, but 
a way out of the difficulty has been found that would improve Wilson's result 
for the Wien effect on equivalent conductance. As a preparation for this line 
of investigation the present author has described a tutorial review of the so-far 
unpublished theory up to the point where the formal formulas for the axial ve- 
locity and ionic field are presented. The details of derivations are given here for 
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the formulas which are not available in the existing scientific journal literature. 
The formal formulas mentioned will be treated in full by using an alternative 
method of evaluation elsewherepS which should produce exact formulas for the 
electrophoretic and relaxation time effects for the Wien effect on conductance 
of binary electrolytes subjected to high external electric field. The results we 
obtain shows, despite the difficulty in the original evaluation of the integrals by 
Wilson, his electrophoretic effect appears to be the lowest order contribution, 



provided that the contributions 



-eXK/8iri]o {np) — 3eX k / '^irrjQ (up) 



can be 



set equal to zero or a finite value in some limit. What such a condition or 
limit would be is not clear at this point unless the velocity expression is fully 
evaluated in (x, p, £) space. Only then will Wilson's result would be properly 
understood. This aspect will be the subject of the forthcoming work. 
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